Understanding the symmetry energy using data from the ALADIN-2000 
Collaboration taken at the GSI Large Neutron Detector 
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The present study deals with the extraction of the symmetry energy from heavy-ion collisions 
at intermediate energies. Using the isospin quantum molecular dynamical (IQMD) model, the 
dependence of the sum of the charge number for fragments with Z > 2 (Abound) on the multiplicity of 
neutrons (M n ) from the projectile spectator fragmentation of 124 Sn and 124 La at 600 MeV/nucleon 
is compared with the experimental results of the ALADIN-2000 Collaboration. The comparison 
suggest a soft symmetry energy. In addition, the sensitivities of the symmetry energy towards the 
Abound dependence on proton multiplicity (M p ), neutron-to-proton single [R(n/p)] and double ratio 
[Rd (n/p)] are also examined. The Zbound dependence of R(n/p) is found to be the most sensitive 
observable towards the symmetry energy. The ALADIN Collaboration should extend the results for 
R(n/p) in the near future. 

PACS numbers: 21.65.Ef, 21.65. Cd, 25.70.Pq, 25.70.-z 



There is an ongoing effort, within the nuclear physics 
and nuclear astrophysics communities, to constrain the 
density dependence of the symmetry energy (DDSE) [1- 
3]. Such constraints are also urgently needed because 
of the unique advantage of the DDSE in understand- 
ing the nuclear equation of state (NEOS) of asymmet- 
ric nuclear matter. In the last couple of years [2-13], 
heavy-ion collisions (HICs) at intermediate energies, due 
to their ability to cover the low-density as well as high- 
density region, have been considered an important tool 
for the determination of the sub- and supra-saturation 
density dependence of the symmetry energy. At present, 
the high-density behavior of the symmetry energy is of 
considerably greater uncertainty compared to the low- 
density behavior. However, a detailed analysis is still 
needed at low density in order to determine the specific 
characteristics such as the nucleon-nucleon cross section, 
the symmetry energy coefficient, method of clusteriza- 
tion etc. [4-6, 14]. Further progress in this direction 
will depend on the results from experiments that are still 
ongoing at laboratories worldwide. 

At present, the results for the high-density behavior 
of the symmetry energy are interesting but vary drasti- 
cally between models. In the following, we will give a 
brief overview of the theoretical as well as the combined 
theoretical and experimental progress. Theoretically, the 
single and double neutron-to-proton ratios [15, 16], sin- 
gle and double tt~ /tt + ratios [7-9, 15, 17], ratio 
[17], K~ I ' K + ratio [7], and isospin fractionation [16] are 
available as sensitive probes to study the high-density be- 
havior of the symmetry energy. In addition, flow param- 
eters are relevant. Flow in heavy-ion collisions is studied 
in terms of the relative and differential directed flow of 
neutron-to-proton, 3 H to 3 He [15], ratio and the differ- 
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ence [10, 11] of neutron-to-proton elliptic flow. Recently, 
the transverse emission of isospin ratios has also been 
suggested as a potential candidate for the determination 
of the symmetry energy [12]. 

In spite of considerable efforts, it is necessary to 
compare the theory with experimental findings. For 
example, in 2009, results from the FOPI collabora- 
tion for the 7r~ /tt + ratio were compared with isospin- 
dependent Boltzmann-Uehling-Uhlenbeck (IBUU04) [8] 
and Lanzhou quantum molecular dynamics (LQMD) [9] 
model calculations. The results were in contradiction be- 
tween the models. In 2011, neutron and proton elliptic 
flow at 400 MeV/nucleon was compared with findings 
from the FOPI Collaboration by two different groups 
[10, 11]. One group predicted the softness of the sym- 
metry energy with 7j = 0.9 using the ratio of neutron-to- 
proton elliptic flow [10]. The other group [11] also found 
evidence for a soft symmetry energy with x = 2 using the 
difference of ncutron-to-proton elliptic flow. However, 
the approaches for the determination of the symmetry 
energy were different: in the first study, the symmetry 
energy was momentum independent, while in the latter 
one, it was obtained from momentum-dependent interac- 
tions. At the end of 2011 [13], the multiplicity of neu- 
trons (M„) data from projectile spectator fragmentation 
at 600 MeV/nucleon at GSI was compared with statisti- 
cal multifragmentation model (SMM) calculations. The 
study pointed out the essential reduction in the strength 
of the symmetry energy. Even after the reduction of the 
symmetry energy coefficient, the data was reproduced 
satisfactorily only in the range Zbound/Zprojectiie > 0.5. 

It is clear that the statistical model failed to reproduce 
the data at low Abound- This raises the important ques- 
tion whether a dynamical model provides an adequate 
description because the mechanism behind the statisti- 
cal and dynamical models is different. In the present pa- 
per, we will try to answer this question by comparing the 
above experimental findings with isospin quantum molec- 
ular dynamical (IQMD) model [18] calculations. In the 
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literature, no one has compared this data with dynami- 
cal model calculations in order to extract the symmetry 
energy. Success in this regard, to a certain extent, may 
help the community reduce the uncertainty in the deter- 
mination of the symmetry energy. In addition, we will 
also check the sensitivity of the Zbound dependence on 
the proton multiplicity(M p ), and the neutron-to-proton 
single [R(n/p)] and double ratios [Roin/p)] with regard 
to the symmetry energy. This will help in providing the 
most sensitive observable for the symmetry energy. 

The study is performed within the framework of IQMD 
model. It has been applied successfully for the determi- 
nation of isospin asymmetric nuclear matter equation of 
state. The isospin effects in the model have been included 
in terms of the symmetry potential, Coulomb potential, 
and isospin dependent nucleon-nucleon (NN) cross sec- 
tions. In the model [18], nucleons are represented by 
the wave packets [19]. The time evolution of the nucle- 
ons, which are under the self-consistency generated mean 
field, is described by the well known Hamilton's equations 
of motion. The Hamilton (H) is the sum of the kinetic 
energy and effective interaction potential. The interac- 
tion potential is composed of Coulomb (Vcoui)i Yukawa 
(Wukawa), local (Vj oc ) and momentum dependent interac- 
tions (Vmdi)- The expressions for Vcoul an d VVukawa have 
been derived by us and others in the Refs. [18]. The local 
interaction potential Moc is originated from the Skyrme 
energy density function. On the basis of this, the local 
potential energy density is expanded as: 



tial in nuclear matter, which is written as follow: 
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where a, /?, and 7 are the parametrized values to specify 
the particular NEOS. The detailed table of the values is 
presented in the Ref. [18, 19]. 

The Eg° m is the symmetry potential energy, which is 
parametrized using the microscopic or phenomenological 
many body theory calculations. The total symmetry en- 
ergy per nucleon employed in the simulation is the sum 
of the kinetic and potential terms and is given as 
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metry kinetic energy coefficient and C. 
is the symmetry potential energy coefficient. 

The momentum dependent potential has been im- 
plemented from Ref. [19] and is expressed as following: 
V M m = C mom ln 2 [e(Ap) 2 + l]j- Q 8{r' - r). Here C mom = 

1.57 MeV and e = 5 x 10~ 4 ,/v 2 ■ The momentum is 
given in units of McV/c. 

Finally, combining all the potentials, we got an isospin, 
density and momentum dependent single particle potcn- 
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In the present simulations, the parameters a, /3, and 7 
are -390 MeV, 320 MeV and 1.14, respectively. They give 
rise to the soft momentum dependent equation of state 
(SMD). The ji parameter in Eq. 2 describes the strength 
of the symmetry energy. The 7; = 0.5 and 1.5 refe 
to the soft and hard symmetry energies, resepectively. 
Moreover, the isospin and energy dependent NN cross 
section in the collision term and the quantum feature in 
terms of Pauli blocking is implemented. 

Following the experimental technique [13], the neutron 
and proton multiplicity is calculated using the cluster 
recognization method. In this method, the particles with 
relative momentum smaller than i"Fcrmi and relative dis- 
tance smaller than Rq are coalesced into one cluster. In 
the present work, the value of Rq and .FWmi are 3.5 fm 
and 268 MeV/c, respectively. 

In order to perform the analysis, we have simulated 
thousands of event for the reactions of 124 Sn+ na 'Sn and 
124 La+ na 'Sn at incident energy 600 MeV/nucleon along 
the whole Zb OU nd range. The effect of different targets 
on the projectile fragmentation is also observed by sim- 
ulating the reaction of 124 Sn+ 124 Sn. In the figures, the 
horizontal (vertical) bar around the displayed symbols is 
the actual Abound bin size (error bar as 95% confidence 
level). The 95% confidence level is ± Standard Devia- 
tion (S.D.). The different colored lines with the same 
colored symbols are drawn spline to guide the trend of 
the different data points as well as reader's eyes. 

It is important to realize the following experimental 
filters in the theoretical simulations: (1) the rapidity dis- 
tribution cut for the projectile spectator fragmentation 
region, and (2) the correlation between the Abound and 
impact parameter (b) [20]. For the first factor, in Fig. 
1(a), the Abound dependence of the rapidity distribution 
of neutrons for 124 Sn+ nat Sn at 600 MeV/nucleon is dis- 
played. At low Abound, a single Gaussian shape distribu- 
tion is observed, which is found to split into two Gaus- 
sians at higher Abound- From here and with the help of 
literature [21], the mid-rapidity region is specified around 
the zero value of rapidity. It is changing into projectile 
(target)-like matter towards the negative (positive) ra- 
pidity values. Using the experimental technique [13], the 
projectile spectator region in the theoretical calculations 
is found to be around -0.9 rapidity value. For the second 
factor, in Fig. 1(b), the impact parameter dependence of 
Abound is displayed. The Z^ oun d is calculated at different 
absolute impact parameters. The Zb oun d increases lin- 
early with the impact parameter. It illustrates that the 
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FIG. 1: (Color online) (a) Rapidity distribution of neutrons 
for 124 Sn+ na 'Sn at 600 MeV/nucleon. The different lines and 
symbols correspond the values at different Zb oun d- (b) Impact 
parameter dependence of Abound with soft (solid circle with 
black line) and hard (open square with red line) symmetry 
energy. The inset is the impact parameter dependence of the 
width of Abound distribution. 



correlation between the two parameters in the dynami- 
cal model agrees with the experimental viewpoint. After 
the above confirmation, in the rest of the analysis, the 
■Zbound is divided into different bins with the size ±2.5 
(Fig. 2-5) by considering the triangular distribution of 
the impact parameter. 

In Fig. 1(b), the sensitivity of the symmetry energy 
on the impact parameter dependence of Abound is also 
checked. No sensitivity is observed at peripheral geome- 
tries, while the visible sensitivity has been found below 
semi-central geometries. This provides an indication of 
the larger Abound with the hard symmetry energy. It 
is possible because of different mechanism for Abound at 
central and peripheral geometries [22]. To assure that 
the visible sensitivity of the symmetry energy is whether 
or not due to the statistical error, the concept of width 
of Abound distribution is used. The width of distribution 
is full width at half maximal (FWHM), which is calcu- 
lated as 2.35x S.D.. The inset in the figure displays the 
impact parameter dependence of width of Abound distri- 
bution with the soft and hard symmetry energy. The 
symmetry energy sensitivity of width of Abound distribu- 
tion resembles with the Zbound distribution, indicating 
sensitivity of impact parameter selector (Zbound) to the 
symmetry energy is due to the statistical fluctuations. 
This is understandable as Zbound is a global variable, 
where the symmetry energy sensitivity becomes weak. 
In the following, we need to search for some potential 
sensitive observables. 

In the left panels of Fig. 2, the Abound dependence 
of M n from the projectile spectator decay is displayed. 
The first two panels in the figure are compared with 
the experimental data of the ALADIN-2000 Collabo- 
ration [13]. Identical to the experimental data, the 
neutron-rich projectile exhibits a maximum of neutrons 
M„ re 11 for 124 Sn and M n re 9 for 124 La at the value 
of Zbound /^Projectile ~ 0.7 in the presence of the soft 
symmetry energy. Above this Abound, more neutrons are 
bound inside fragments and hence the multiplicity of neu- 
trons decreases. 
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FIG. 2: (Color online) Abound dependence of M n (M p ) in left 
(right) panels from projectile spectator fragmentation for dif- 
ferent projectile-target combinations. The open square (solid 
circle) including lines have the same meaning as in the (b) 
part of the Fig. 1. The experimental data of ALADIN Col- 
laborations is displayed by blue stars. 



Interestingly, the multiplicity of neutrons presents sen- 
sitivity to the strength of the symmetry energy. In all the 
panels, the soft symmetry energy yields more neutrons 
compared to the hard symmetry energy. The sensitivity 
towards the symmetry energy has been seen more clearly 
from the more neutron-rich projectile 124 Sn compared to 
the neutron-poor 124 La. The previous comparison with 
the SMM calculations reproduced the data satisfactorily 
only in the range Zbound/ZprojectUe > 0.5. The present 
IQMD model reproduces the experimental data along the 
whole Zbound range. However, some discrepancies are 
observed for neutron-poor projectile 12 La compared to 
neutron-rich projectile 124 Sn. This depends on the dif- 
ferent strength of the symmetry energy and Coulomb in- 
teractions in different projectiles. Owing to the more 
neutron- rich nature of 124 Sn, stronger symmetry energy 
effect has been observed and hence a good agreement 
with the data is obtained. While in case of 124 La, due to 
the presence of 4 extra protons compared to 124 Sn, the 
contribution from Coulomb interactions results in the de- 
viation of the calculations from the experimental data. 
The enhanced production of protons in Fig. 2(d) com- 
pared to Fig. 2(b) is further supporting the discrepancies. 

In addition, the multiplicity of neutrons from the pro- 
jectile spectator fragmentation is strongly dependent on 
the type of projectile [(a) and (c)], while, is weakly depen- 
dent on the type of the target [(a) and (e)]. The related 
study has also been conducted at relatively low energy 25 
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FIG. 3: (Color online) Abound dependence of R(n/p) for 
neutron-rich (neutron-poor) 124 Sn ( 124 La) projectile specta- 
tor fragmentation in top (bottom) panels. The symbol and 
lines have the same meaning as in Fig. 2. The vertical dashed 
lines correspond to rough peak positions of Zb oun d for R(n/p) 
with hard and soft symmetry energy, respectively (from left 
to right). 



MeV/nucleon for 40 Ca + 40 < 48 Ca, 46 Ti reactions at the 
4?r CHIMERA detector [23] . They showed that competi- 
tion between fusion-like and binary reactions in selected 
ccntrality bin can constrain the parametrization of the 
symmetry energy. 

In the right panels of Fig. 2, the Abound dependence 
of proton multiplicity M p from the projectile spectator 
decay is displayed. The behavior of the curves resembles 
with the left panels. The differences between the left and 
right panels are as follows: (1) The M p is more with the 
hard symmetry energy, which was true with the soft for 
the M n . This is because of the repulsive (attractive) na- 
ture of symmetry energy for neutrons (protons). At the 
freeze-out time, the soft symmetry energy has high mag- 
nitude, which will lead to the production of neutrons. It 
has alternatively led to the less production of protons; 

(2) More protons are produced for neutron-poor projec- 
tile 124 La compared to the neutron-rich projectile 124 Sn; 

(3) The sensitivity of symmetry energy is weak on the 
M p in 124 La projectile, which is due to the dominance 
of Coulomb interactions. These factors are revealing the 
importance of the symmetry energy for neutron-rich re- 
actions. 

Now we will analyze the Abound dependence of neutron- 
to-proton single ratio R(n/p) from the projectile spec- 
tator decay. In Fig. 3, the R(n/p) is displayed for 
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FIG. 4: (Color online) Zbound dependence of double neutrons 
to protons ratio (R D (n/p) = l °J r i^f" ) with soft and 

hard symmetry energy. 



the 124 Sn+ nat Sn (top panel) and 124 La+" at Sn (bottom 
panel) reactions. The single ratio is found to be more sen- 
sitive towards the soft symmetry energy. Particularly, the 
R(n/p) with the soft (hard) symmetry energy is higher 
(lower) than the ratio of system [shown by blue horizontal 
line] . This is attributable to the large (small) magnitude 
of soft (hard) symmetry energy at freeze-out time. In ad- 
dition, Figs. 2 and 3 has some similarities and differences. 
The behavior of the single ratio curves almost resemble 
with the Zbound dependence of M n and M p , i.e. it also 
exhibits a maximum or saturates with Abound- The dif- 
ference is that the maximal value for R(n/p) is noticed 
at different Zbound/Zprojectiic with the soft (at 0.7) and 
hard (at 0.5) symmetry energy, which was at the same 
Zhound/Z pro j ec tiie for the maximal M n and M p (at « 0.7). 
It means that the single ratio has more (less) neutron- 
rich effects along Abound with the soft (hard) symmetry 
energy. This is once again due to different strength of 
symmetry energy and Coulomb interactions originated 
from the soft and stiff parametrization of the symmetry 
energy. 

As discussed in the Ref. [4], the double neutron-to- 
proton ratio Roin/p) from the two isotopes of a system 
is the most sensitive observable for the determination of 
the symmetry energy. In the present study, the projec- 
tiles are isobars and not the isotopes. The Abound depen- 
dence of Rr>{n/p) using the isobars is shown in Fig. 4. 
From the comparison of Fig. 3 with Fig. 4, it is found 
that the Abound dependence of the single as well as dou- 
ble ratio is more sensitive to the soft symmetry energy. 
In contrary, the double ratio, comparing to the single 
ratio, has relatively weak sensitivity toward the symme- 
try energy along the whole inbound region. Moreover, the 
Ru{n/p) with the hard symmetry energy is weakly de- 
pendent on the Abound- This may be the effect of four 
extra protons in the 124 La projectile, where the contribu- 
tion from the Coulomb interactions can not be canceled, 
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FIG. 5: (Color online) Ztound dependence of sensitivity of 
different observables (shown in Fig. 2, 3, and 4) towards 
the symmetry energy. (a) for 124 Sn+ nai Sn and (b) for 
124 La+ nQ *Sn. 



even by taking the double ratio. 

Finally, to declare the most sensitive observable to- 
ward the symmetry energy among the all three, i.e. M n , 
R(n/p), and Rd(u/p), we have defined the sensitivity 
factor (St) as follow 



St= 1^=0.5 -^7«=1.5| xlQ0 



(4) 



where X = M n , R(n/p) or R£>(n/p). 

From the results of St (Fig. 5), we can see that the 
St values for all the parameters are more sensitive for 
neutron-rich projectile, illustrating that neutron-rich nu- 



cleus is a good candidate for the extraction of the symme- 
try energy. As was expected from Fig. 4, the St factor 
of Rn(n/p) shows the least sensitivity to the symme- 
try energy in this particular situation of different iso- 
bars. Between the other two, i.e. M n and R(n/p), the 
R(n/p) shows more sensitivity to the symmetry energy. 
This sensitivity becomes more stronger with the increase 
of Abound- Therefore, it is advisable that experimental 
Collaboration should also make the data for the R(n/p) 
available, which can help to check the validity of our re- 
sults for R(n/p) and Rjj(n/p). 

In conclusion, using the IQMD model, the ALADIN- 
2000 Collaboration data for multiplicity of neutrons ex- 
tracted from the projectile spectator fragmentation is 
satisfactorily reproduced with the soft symmetry energy. 
Along with multiplicity of neutrons, proton multiplicity 
has also a peak around ■Zbound/^porjectiic ~ 0-7 with the 
soft as well as the hard symmetry energy. Interestingly, 
the R(n/p) produced a peak around ^bound/^porjectiie ~ 
0.7 (0.5) with the soft (hard) symmetry energy. The 
R(n/p) for the different isobar projectiles is the most sen- 
sitive observable towards the symmetry energy in com- 
parison to multiplicity as well as Ro(n/p). In the near 
future, we hope experimental Collaborations make the 
results for proton multiplicity and R(n/p) available. 
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